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The capability of the OVERFLOW code to accurately compute high-speed turbulent 
boundary layers and turbulent shock-boundary layer interactions is being evaluated. Con- 
figurations being investigated include a Mach 2.87 flat plate to compare experimental ve- 
locity profiles and boundary layer growth, a Mach 6 flat plate to compare experimental 
surface heat transfer, a direct numerical simulation (DNS) at Mach 2.25 for turbulent quan- 
tities, and several Mach 3 compression ramps to compare computations of shock-boundary 
layer interactions to experimental laser doppler velocimetry (LDV) data and hot-wire data. 
The present paper describes outlines the study and presents preliminary results for two of 
the flat plate cases and two small-angle compression corner test cases. 


Nomenclature 


Cf Local skin friction coefficient 

M Local Mach number 

V Velocity 

S Boundary layer thickness 

S* Boundary layer displacement thickness 

0 Boundary layer momentum thickness 

P s Static pressure 

T Temperature 

x Streamwise coordinate 

y Crossflow coordinate 

z Wall-normal coordinate 

L sep Separation length 

Subscripts 

R e f Reference value for data reduction 

o Stagnation condition 


I. Introduction 

A useful thermal analysis tool for the space shuttle and future space vehicles needs to accurately model 
turbulent boundary layers in order to predict the heat transfer on the surface. Shock-boundary layer 
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interactions cause a significant increase in wall heating in the reattachment region of a shock-separated 
boundary layer, so they must be accurately modeled as well. The present work is carried out to continue 
the development of the capability to model high-speed compressible flows over complex vehicles. The goal 
of the current project is to evaluate the performance of the turbulence models in the OVERFLOW code 1,2 
for high-speed, non-reacting flow. 

The OVERFLOW code is an overset (chimera) grid Navier-Stokes solver, which makes it ideal for com- 
puting the flow fields over complex geometries. 2 3 The code utilizes the Swanson/Turkel matrix dissipation 
model 4 which, with the correct parameters, 5 allows the central- differenced scheme to mimic a Total Variation 
Diminishing (TVD) upwind biased scheme in regions of flow discontinuities. This method takes advantage 
of the speed of central differencing while maintaining the accuracy of the upwind TVD scheme. The idea 
of using OVERFLOW as a thermal analysis tool has received some criticism based on the fact that it was 
designed to compute aerodynamic forces in transonic and low supersonic flows. Previous work 5,6 has indi- 
cated that it is indeed capable of accurately modeling laminar hypersonic flows and does an adequate job of 
capturing heat transfer at those speeds. The present study is now directing attention at validating the use 
of OVERFLOW for high-speed turbulent flows. 

Three turbulence models in OVERFLOW will be considered in this paper: the Spalart-Allmaras model, 7 
the SST model, 8 and the Lag model. 9 Results from each model will be compared to several different 
experiments and direct numerical simulation (DNS) computations in order to benchmark their behavior. 

A. Code Description 

OVERFLOW is a structured (overset) grid Navier-Stokes flow solver. It uses a finite- difference formulation 
with flow quantities stored at the grid nodes. OVERFLOW has central- and Roe upwind- difference options, 
and uses a diagonalized, implicit approximate factorization scheme for time advancement. Local timestep- 
ping, multigrid techniques and grid sequencing are all used to accelerate convergence to a steady state 
solution. In this study, a 2nd-order central differencing with Swanson/Turkel type 2nd/4th-order matrix 
dissipation is utilized. For a complete discussion of the scheme, see Olsen et al. 5 

B. Turbulence Models 

OVERFLOW has numerous turbulence models coded in, but this study will focus on three models, the 
Spalart-Allmaras (SA) model, 7 the SST model, 8 and the Lag model. 9 The turbulence models included in 
the code all make use of the Boussinesq approximation, meaning that the effects of the Reynolds stress 
terms are included in the Navier-Stokes equations through an eddy viscosity. The SA model is a popular 
model for aerodynamic flows that solves a single partial differential equation directly for the eddy viscosity. 
The SST model is a two-equation model which is a hybrid of the k-e and k-o; turbulence models and solves 
equations for the turbulent kinetic energy and the specific dissipation rate. The eddy viscosity is then derived 
from an algebraic relation that chooses the minimum value between the standard formulation based on the 
turbulence kinetic energy and the specific dissipation rate and a formulation based on an assumed value of 
the normalized Reynolds shear stress. The Lag model is a recent extension of the SST model that accounts 
for non-equilibrium effects by carrying an additional differential equation that relaxes the eddy viscosity to 
the equilibrium value. This third equation accounts for the time required for the turbulence to respond to 
changes in the mean flowfield. Each of the test cases is computed with all three turbulence models to allow 
comparisons of the models. 

C. Study Objectives 

The objective of the current study is to benchmark the performance of these turbulence models, as imple- 
mented within OVERFLOW, for the high-supersonic to hypersonic regimes experienced by the space shuttle 
during ascent. Surface properties, such as surface pressure and heat transfer, are the primary variables of 
interest, however to aid in the development and implementation of improved models and model corrections, 
profiles of flow variables and turbulent quantities are also being generated. 

This paper presents the surface properties and select profiles for a few of the test cases that have been 
completed. There is much work that remains to be done in this study, however in the interest of generating 
useful feedback, the entire project will be outlined. 
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II. Test Cases 


The test cases are each chosen to address a primary objective in this study. In order to validate the 
code’s ability to capture a zero pressure gradient boundary layer, several flat plate cases will be run. Two- 
dimensional compression ramp experiments have been selected to analyze the shock/boundary layer inter- 
action behavior, namely the size of the separation region, the profiles of turbulent quantities downstream of 
the shock, and the quality of the prediction in the recovering boundary layer. Since experiments often have 
to make some compromises in terms of what measurements can be made, multiple experiments have been 
chosen for both geometries to capture specific features of the turbulent boundary layers. The experiments 
have been chosen based on perceived quality of the data and Mach number of interest. 

A. Flat Plate 

In order to adequately address each relevant feature of a high-speed boundary layer, four datasets have been 
selected for comparison. An experimental study by Cary 10 at Mach 6 has been chosen for its surface heating 
data. A DNS simulation by Gatski 11 provides two detailed profiles of mean and turbulent properties. Finally, 
an experiment conducted by Smits, et al. 12 at Mach 2.87 has profiles at multiple locations upstream of an 8° 
compression corner, and will be used for comparison of high Reynolds number experimental velocity profiles 
and boundary layer growth measurements. 

B. Compression Ramp 

Two supersonic compression ramp datasets have been chosen for comparison. The Mach 2.9 experiments of 
Kuntz 13, 14 will provide two component laser doppler velocimetry (LDV) measurements of turbulent boundary 
layers on ramps of 8, 12, 16, 20, and 24 degrees, and the Mach 2.9 series of experiments 12, 15 conducted at the 
Princeton Gas Dynamics Laboratory will be used to compare turbulence quantities downstream of the shock 
and for skin friction and surface pressure distributions on ramps of 8, 16, 20, and 24 degrees. Settles and 
Dodson reported in their 1994 shock-boundary layer interaction database 16 that these two datasets disagree 
on the magnitude of the Reynolds stresses by as much as a factor of 4. To date, the authors are not aware of 
an explanation of this discrepancy, so following the suggestion of Settles and Dodson, these values are taken 
to be the bounds of the actual Reynolds stresses. 

The Mach 3.0 LES simulations of Rizzetta et al. 17,18 are also being considered for comparison, however 
difficulties have been encountered while trying to compare to this dataset and it is unknown if these issues 
will be worked out. 


III. Solution Methods 

A. Grid Generation and Boundary Conditions 

The flat plate test cases use simple single-zone grids that do not model the radius of the plate leading edge, 
and the viscous wall extends all the way to the inflow plane. The compression ramp cases use a system 
similar to that depicted in Figure 1. The inflow turbulent velocity profiles given in each experimental case is 
generated on the long flat plate region upstream of the corner, which is computed together with the ramp as 
one computation. The OVERFLOW cases are run fully turbulent, so the length of this plate must changed 
to specify the ’transition point’. The length of the flat plate is chosen such that the momentum thickness 
in matched to the experiment as closely as possible. The corner and ramp regions are defined by three 
grids: two high resolution near-wall grids (split for load-balancing on multiple processors) and a grid in the 
freestream. It was found that the hyperbolic tangent stretching function used for the wall-normal direction 
did not place the desired number of grid points in the freestream above the shock wave, so the freestream 
was given its own grid. The final grid is a low resolution sponge grid to contain the leading edge shock wave 
and prevent it from reflecting back and interfering with the boundary layer development. 

The boundary conditions for both grid systems are freestream at the inflow boundary, characteristic 
freestream on the boundary opposite the viscous wall, and extrapolation at the outflow. The viscous wall is 
either adiabatic or isothermal, based on the specific dataset being compared to. 

It was discovered that adding the flat-top to the compression ramp downstream of the experimental data 
changed the size of the separation bubble for a large ramp angle. This was taken to mean that disturbances 
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from the 0 th order extrapolation outflow boundary condition were propagating upstream through the subsonic 
portion of the boundary layer. While it has not yet been determined if this holds true for smaller ramp angles, 
the flat tops were also added to these cases for the sake of consistency. 
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Figure 1. Grid topology for 24° ramp with every other grid line shown, (a) Overall grid system, (b) Detail 
view of corner region. 


B. Solution Method 

Each of the solutions is initialized to freestream conditions and iterated on a sequence of progressively finer 
grids until the finest grid level is obtained. Coarse grids are created internally in OVERFLOW by specifying 
a number of grid levels (three grid levels would be the original fine grid, a medium grid with every other 
grid point used, and a coarse grid with every fourth grid point used). For this study, three grid levels were 
used and between 500 and 1,000 iterations were performed on the coarse and medium grid levels to startup 
the solution. 
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In order to rapidly obtain a high quality steady-state solution, the dissipation parameters for the matrix 
dissipation scheme 4 are ramped down to their optimum level as the solution converges. The four parameters 
controlled are the 2nd and 4th-order dissipation coefficients, and ^4, and the eigenvalue limiters V Cr? and 
V e£ (DIS2, DIS4, VEPSN, and VEPSL in OVERFLOW, respectively). The presented solutions were started 
with these values set to (10, 0.2, 1.0, 1.0) and were ramped down to (2, 0.1, 0.3, 0.3). For more information 
on this process, see Lillard et al. 5 6 

The solutions are considered to be converged with the optimal dissipation settings when: (1) the residuals 
have dropped by more than three orders of magnitude and, (2) the skin friction at several representative 
locations has converged to a constant value. Grid convergence was determined using the Richardson ex- 
trapolation similar to the method of Roache. 19 The solutions in the present paper were considered to be 
adequately grid resolved when the average difference between the high resolution skin friction and the ex- 
trapolated skin friction in the corner /ramp region was less than 3%. Grid convergence results are presented 
in Table 1 and Table 3. 

Table 1. Inflow matching and grid convergence results for compression ramp test cases 


Case 

e 

@ Target 

< 5 * 

^ T arqet 

RMS A C f 

Max A Cf 

A L/Sep 

8° Ramp (SST) 

1.001 

1.002 

0.78% 

2.0% 

2.1% 

8° Ramp (SA) 

0.986 

0.993 

1.8% 

8.5% 

7.9% 

8° Ramp (Lag) 

0.994 

0.996 

1.76% 

7.4% 

3.0% 

16° Ramp (SST) 

0.999 

1.060 

2.2% 

17% 

2.6% 

16° Ramp (SA) 

0.995 

1.063 

2.1% 

8.8% 

0.58% 

16° Ramp (Lag) 

1.015 

1.078 

2.9% 

14.4% 

0.77% 


IV. Preliminary Results 

Work on this project is continuing; however, several computations have been completed and are presented 
below. These solutions include the Gatski flat plate and the mean flow portion of the 8° and 16° compression 
corners from the Princeton series of experiments. The reference values in Table 2 are used to scale the output 
of the following data. 


Table 2. Reference values used in data reduction 


Case 

Flat Plate 

8° Ramp 

16° Ramp 

MR e f 

2.25 

2.87 

2.85 

VrcJ (m/s) 

260.9 

590 

580 

S Ref (mm) 

3.8 

26 

26 

PRef (kPa) 

- 

22.52 

23.56 


A. Flat Plate 

Two of the flat plate test cases have been run to date, namely the Gatski DNS test case and the feeder 
plate of the 8° compression ramp experiment from the Princeton series. The data available from the Gatski 
dataset includes skin friction and velocity profiles at two points on the flat plate. Table 3 shows the skin 
friction results for each turbulence model against the DNS value, and Figure 2 shows the two velocity profiles. 

The analysis of the flat plate boundary layer upstream of the 8° compression ramp is given in Table 
tabiccfp. The boundary layer displacement thickness and momentum thickness comparisons seem to be 
within the experimental scatter for all turbulence models (keeping in mind that 0 at x = — 25.4mm was used 
to define the length of the flat plate). The skin friction predictions range from 9% to 13% too large for the 
OVERFLOW solutions. It is difficult to call one turbulence model better than the others based on this data 
alone. 
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Table 3. Gatski flat plate grid convergence and skin friction results 


Case 

RMS A C f 

Max A Cf 

Cf x = 7.8m 

Cf x — 8.8m 

DNS 

- 

- 

2.59 x 10“ 3 

2.25 x KT 3 

SST 

0.53% 

0.55% 

2.01 x 10“ 3 (22%) 

1.97 x KT 3 (12%) 

SA 

0.63% 

0.65% 

2.07 x 10 -3 (20%) 

2.03 x 10 -3 (9.7%) 

Lag 

0.60% 

1.07% 

2.02 x 10“ 3 (22%) 

2.01 x KT 3 (11%) 



v/v Ref 


Figure 2. Velocity profile at x = 8.8in for Gatski flat plate test case. 


B. Compression Ramp 

Freestream conditions for the compression ramp cases shown are given in Table 5. Figure 3 shows the skin 
friction distributions on the 8° and 16° ramps. For both of these cases, the skin friction from the upstream 
boundary layer is too high, and all the OVERFLOW computations indicate separation when it was not 
observed in the experiments. This aside, the Lag model seems to perform better with increasing shock 
strength, while the SA performance decreases. The SST model seems prone to separating too early. 

Figure 4 shows the static pressure in the flowfield and on the surface of the 8° ramp. For this case, all 
of the models do a descent job representing the data. Figure 5, which shows the pressure for the 16° ramp, 
is a little more interesting. It appears that, as expected from Figure 3 above, the SST model separates too 
early and the oversized separation bubble throws off the pressure in the recovering boundary layer. The Lag 
model predicts the location of the pressure rise very well, but the pressure in the recovering boundary layer 
falls a little below the measured pressure. The SA model, on the other hand, misses the pressure rise but 
gets the recovery pressure profile. The reversal in which model better predicted the surface pressure and 
skin friction seems to indicate that these differences are within the experimental uncertainty. 

Looking at the contour plots in Figures 4 and 5, it appears that there are large disturbances and a fat 
shock in the inviscid flowfield behind the shock wave. This is an artifact of the numerical dissipation scheme 
used in OVERFLOW, which is known to smear shocks over 5-7 grid points. To see whether or not this 
was affecting the surface properties, an case was run with a special grid to capture and contain the shock 
wave. This grid, shown in Figure 6a, is aligned closely with the shock predicted from the 16° computation. 
Figure 6b shows the skin friction for the 16° ramp with and without the shock grid. As you can see, the 
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Figure 3. Profiles of skin friction on (a) 8° ramp, (b) 16° ramp. 
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Figure 4. 
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and (b) surface pressure on 8° ramp. 
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x ^Ref 

Figure 5. (a) Contours of and (b) surface pressure on 16° ramp. 
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Table 4. Boundary layer development on feeder flat plate for 8° compression corner case 


Case 

Cf 

6 (mm) 

5* (mm) 

6 (mm) 

x = — 50.8mm 

Experiment 

1.00 x 10 3 

26.0 

6.7 

1.3 

SST 

1.09 x 10 3 

24.5 

6.60 

1.29 

SA 

1.13 x 10 3 

27.1 

6.58 

1.27 

Lag 

1.11 x 10 3 

26.1 

6.61 

1.28 

x = — 25.4mm 

Experiment 

1.00 x 10 3 

26.5 

6.6 

1.3 

SST 

1.09 x 10 3 

24.8 

6.71 

1.30 

SA 

1.13 x 10 3 

27.4 

6.66 

1.28 

Lag 

1.11 x 10 3 

26.4 

6.61 

1.29 

x = — 15.2mm 

Experiment 

0.990 x 10 3 

26.5 

7.3 

1.2 

SST 

1.09 x 10 3 

24.9 

6.75 

1.31 

SA 

1.13 x 10 3 

27.5 

6.69 

1.29 

Lag 

1.11 x 10 3 

26.6 

6.71 

1.30 


Table 5. Freestream and inflow conditions used in compression ramp computations 


Ramp Angle 

8° 

16° 

24° 

Mo© 

2.87 

2.85 

2.84 

To (k) 

280 

268 

262 

Po (kPa) 

680 

690 

690 

S (mm) 

26 

26 

23 

S* (mm) 

6.7 

6.3 

6.1 

S (mm) 

1.3 

1.3 

1.2 

% upstream (mm) 

-25.4 

-25.4 

-63.5 


differences are negligible (less than 1% along the entire length), so the shock grids were left off to ease the 
load of obtaining grid-converged solutions for the present paper. Shock grids will be used once profiles of 
mean and turbulent properties are analyzed in detail. 

Figures 7 and 8 show velocity and local Mach number profiles for the 8° and 16° ramps respectively. For 
the 8° ramp, the shock-boundary layer interaction is weak enough that most of the differences between the 
three turbulence models is small. However, for the 16° ramp, characteristic features of RANS predictions 
begin to appear. In each of the profiles downstream of the shock, the flow near the wall slows down too 
much. The Lag model does a much better job than the other turbulence models in this regard and fits the 
near- wall data very well for profile 5, well within the experimental uncertainty. 

As the present data shows, these two ramp angles do not really create a strong enough interaction to 
significantly separate these three turbulence models; solutions for steeper ramp angles will be required. 
Work is presently underway to obtain steady solutions of the 20° and 24° compression corners, however 
initial predictions have been very poor and often do not converge to a steady solution. One option currently 
being explored is running OVERFLOW in a time-accurate unsteady calculation and taking time-averages 
of the surface properties. Figure 9 shows initial results for skin friction and surface pressure the 24° ramp. 
The SA model converged to a steady solution, however the correlation with the experimental data is very 
poor-worse than would have been expected beforehand. The SST solution shown is the time average of an 
initial unsteady computation. This method looks promising and work is currently underway to make sure 
that the time resolution is sufficiently fine and the run time is long enough to give accurate statistics. 
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Figure 6. 16° ramp with shock grid, (a) Grid system, (b) Contours of p Ps , (c) Computed skin friction with 

^ Ref 

and without shock grid. 
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Figure 7. 8° ramp: (a) Mach contours with profile locations indicated, (b) Mach profiles, (c) Velocity profiles. 
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SA Experiment 


Figure 8. 16° ramp: (a) Mach contours with profile locations indicated, (b) Mach profiles, (c) Velocity profiles. 
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Figure 9. (a) Skin friction and (b) Surface pressures on 24° ramp. Note that SA curves are from a steady-state 
solution and the SST curves are a time averaged unsteady solution. 
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V. Conclusions and Future Work 


The primary objective of this paper is to evaluate the capability of the OVERFLOW code to accurately 
predict heat transfer in turbulent hypersonic flows and identify what aspects of the code need improvement 
to make OVERFLOW a useful aerothermal analysis tool. Multiple flat plate and compression corner ex- 
periments have been selected for the evaluation of specific features of turbulent high-speed boundary layers. 
Preliminary results seem to show that for weak shock-boundary layer interactions, the SA, SST, and Lag 
turbulence models behave similarly, however differences begin to appear for a 16° ramp at Mach 3. All of 
the models predict separation where the experiments reports none, thought the SST model is the worst at 
over predicting the separation size. The Lag model appears to do the best overall of the three models for the 
few cases shown in this paper. Current and future work in this study is focused on obtaining high-quality 
solutions with rigorous grid resolution studies for each test case. It appears that some of the steeper ramp 
angles may require special treatment (i.e. running unsteady), and these methods will be explored. 
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